Micro-SQUID characteristics 
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We report on the dependence on field and temperature of tfie critical current of micro-SQUIDs: 
SQUIDs with diameters as small as 1 micron, using Dayem bridges as weak links. We model these 
SQUIDs by solving the Ginzburg-Landau equations with appropriate boundary conditions to obtain 
the supercurrent-phase relationships. These solutions show that the phase drops and depression 
of the order parameter produced by supercurrent flow are often distributed throughout the micro- 
SQUID structure, rather than being localized in the bridge area, for typical micro-SQUID geometries 
and coherence lengths. The resultant highly non-sinusoidal current-phase relationships Id'p) lead 
to reduced modulation depths and triangular dependences of the micro-SQUID critical currents on 
applied magnetic flux /c("I'). Our modelling agrees well with our measurements on both Al and Nb 
micro-SQUIDs. 
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INTRODUCTION 



Several different types of specialized Superconducting 
Quantum Interference Devices (SQUIDs) have been de- 
veloped for measuring the magnetic response of small 
samples. One such device is the "micro-SQUID" : a thin 
film DC-SQUID with Dayem-Bridges as Josephson Junc- 
tions In the micro-SQUID, the entire device is fabri- 
cated by electron beam lithography, and the SQUID loop 
itself serves as the flux-input coil. Micro-SQUIDs have 
the advantages of very small pickup areas (about 1/im^), 
and relatively small sensitivity to in-plane applied mag- 
netic fields (since they can be made very thin, less than 
20 nm). Micro-SQUIDs have been used for the observa- 
tion of persistent currents in 2 dimensional electron gas 
rings H], for the study of the mechanisms of magnetiza- 
tion reversal in ferromagnetic particles as small as 3 nm in 
diameter |^,^, and have been integrated into a scanning 
SQUID- AFM with high magnetic field spatial resolution 
However, these SQUIDs have two disadvantages: 
they have hysteretic current-voltage characteristics, and 
relatively low modulation depths in their critical current 
- flux characteristics. In this paper we report detailed 
measurements of micro-SQUID characteristics, and com- 
pare them with model calculations. These calculations 
model the micro-SQUID characteristics well, and could 
represent a valuable tool for optimizing the properties of 
this class of SQUIDs. 



MEASUREMENT 

The basic operating properties of micro-SQUIDs are 
understood: They have a hysteretic V(I) characteristic, 
induced by the propagation of a hot spot. As the current 
is ramped up from zero, the micro-SQUID transits from 
the superconducting to the normal state at a critical cur- 
rent Ic. A voltage step is generated as the normal state 
resistance of the junction appears and the dissipated en- 
ergy heats the entire micro-SQUID loop. When the cur- 
rent is lowered the micro-SQUID stays in the resistive 
state until the current is much smaller than Ic. This ther- 
mal hysteresis excludes the usual current biasing schemes 
used for DC-SQUID readout. Therefore a detection tech- 
nique suitable for hysteretic devices |^] is implemented: 
A computer controlled circuit simultaneously triggers a 
current ramp and a 40MHz quartz clock. As soon as a 
dV/dt pulse of a preset height is detected at the micro- 
SQUID, the clock stops and the current is set to zero. 
The clock reading is transferred to the computer, and the 
cycle begins again. The critical current is proportional to 
the duration of the current ramp. The fastest repetition 
rate is lOkHz, limited by the time needed to settle the 
current. A single wire is sufficient to connect the micro- 
SQUID, since the dV/dt pulse is detected on the current 
biasing lead. Every time the critical current is measured, 
the flux state in the micro-SQUID is sampled, and every 
time the micro-SQUID becomes normally conducting the 
external field penetrates. As the current is reduced (in 
40ns) to zero the micro-SQUID structure becomes pro- 
gressively superconducting again and screening currents 
are set up to quantize the total flux through the micro- 
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SQUID. In the limit of high critical currents different flux 
configurations can be stabilized in the micro-SQUID ring 
during the backswitching, leading to multivalued Ic vs B 
characteristics. 

However, the details of the dependence of the critical 
current of the micro-SQUID on magnetic field and tem- 
perature have not been well characterized. We present 
in this manuscript a detailed description of the underly- 
ing physics necessary for the understanding and improve- 
ment of this type of SQUID. In this study we made Ic(T) 
and Ic vs B measurements on micro-SQUID devices made 
of aluminum or niobium. Each SQUID consists of a 1/im 
square loop. Bulk aluminum and niobium have very dif- 
ferent superconducting properties, such as the coherence 
length ^, and the normal to superconducting transition 
temperature Tc. The characteristics of micro-SQUIDs 
(see e.g. Fig.'s 2 and 4) do not resemble those of ideal 
Josephson SQUIDs, in that they have relatively shallow 
modulation depths and triangular Ic vs B interference 
patterns at low temperatures. 

We believe that these non-ideal characteristics occur 
because the weak links in these SQUIDs are Dayem- 
Bridges, with dimensions comparable to the coherence 
length. It is well known that "long" Josephson weak 
links can have non-sinusoidal Josephson current-phase 
relationships and that non-sinusoidal Josephson 

current-phase relationships can in turn lead to non-ideal 
SQUID Ic vs B interference patterns Faucher et 

al. have attributed the non-ideal critical current vs. 
flux characteristics of their micro-SQUIDs to the large 
kinetic inductance of the micro-bridges. We believe the 
treatment given here is equivalent to that of Faucher et 
al. in terms of the properties of the micro-bridges, but 
it goes further, in that it includes the properties of the 
entire SQUID loop itself. We model the SQUID current- 
phase relationships by solving the 2-d Ginzburg-Landau 
equations. It is important to solve these equations in the 
full geometry of the micro-SQUID, since as previously 
pointed out in general for Dayem-Bridges 0, and as is 
born out by our modelling in this particular case, the 
phase drops and order-parameter depression associated 
with supercurrent flow often extend far from the micro- 
bridge. 



MODEL 



scribing the superconducting state, and A is the vector 
potential. The supercurrent J is given by: 



Amri m*c 



(2) 



For the purposes of the paper, we include the effects of 
the vector potential A in a lumped circuit element model 
(see Eq. |l^). Therefore in our solution of the GL equa- 
tions we set ^ = 0. Writing the complex order parame- 
ter as V' =1 "0 I 6*''') using a coordinate system in which 
the SQUID is in the xy plane, and neglecting the z de- 
pendence of the gradients of tp, the real (Eq. |^) and 
imaginary (Eq. ^) parts of Eq. |l| become: 
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and 
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(4) 



Setting I ■(/) 1= / I ?/) loo, where / is a real function of x 
and y, and | i/; \aa is the unperturbed value of the order 
parameter, /? | i/^ -a = /2m*£,'^{T) pL and using 
reduced units x' — x/^{T),y' = y/£,{T), Eq. becomes 



dx'-^ dx''^ 



Eq. ^ becomes 
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The first Ginzburg-Landau (GL) differential equation The supercurrent density J becomes 
is fill 



aip + P \ ip l"^ ip + 



2m* 



!v-f:ii (1) 



J : 



e*h 



I ^ L f 



m*^{T) 



^ dip ^ dip 
dx' ^ dy' 



(7) 



where m* = 2m and e* = 2e are the mass and charge of It can be shown that Eq. ^ is equivalent to setting the 
the Cooper pair, ip is the complex order parameter de- divergence of the supercurrent (Eq. ^ equal to zero. 
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FIG. 1. Solution to the 2-dimensional Ginzburg-Landau 
differential equations for the quantum mechanical phase Lp and 
normalized superfluid density / for a micro-SQUID. The ge- 
ometry of the SQUID is that of the Al micro-SQUID for which 
measurements are presented in Fig. ^. The length of the 
bridge S — 290nm, bridge width=66 nm, the SQUID length 
(horizontal, x) is 1.3/im and width (vertical, y) is lAbfim. The 
calculations for this figure were done for ^=580nm, with a to- 
tal phase drop across the SQUID of (ps ~ t. Greyscale images 
of the superconducting phase ip (a) and superfluid density / 
(b) in the SQUID are shown. Shown in (d) are cross-sections 
through these images, along the path shown in (c). For this 
value of S/^(r)=0.5, the phase drop and superfluid depres- 
sion are mostly localized in the micro-bridge region. 

By taking discrete steps S in both the x and y direc- 
tions, the differential equation Eq. || can be cast as a 
difference equation 



+ fi-l,j + f'Lj + l + /tj- 



~^,,,f+6^)f,,,^S'fl^, (8) 



(11) 

We solved Eq.'s ^ and ^ using a general non-hnear dif- 
ferential equation solver subroutine (SOR) using the suc- 
cessive over-relaxation method with Chebyshev acceler- 
ation An initial guess was made for the matrices 
fij and ipi,j, SOR was used to search for a solution for 
/i.j of Eq. 1^ for a pre-determined number of iterations, 
with ipij fixed. The result for fij was then fixed, and a 
solution for ipi^ in Eq. ^ was sought for the same num- 
ber of iterations. This procedure was iterated until the 
sum of the absolute values of the deviations from equal- 
ity in Eq.'s |^ and ^ were both less than a fixed error sum 
value. For the results reported here, the error sum value 
was chosen to be 10~^. This procedure gave results in 
agreement with those reported by Likharev and Yakob- 
son for long 1-D microbridges For modelling of our 
micro-SQUIDs we chose the boundary conditions: 1) (p 
was fixed at along the entrance to the micro-SQUID 
structure (at the left of Fig. |l|a)), and fixed at ips at the 
exit of the micro-SQUID structure (to the right); 2) / 
was chosen to be 1 at both the entrance and exits to the 
micro-SQUID; and 3) the components of the gradients 
of / and If were taken to be zero normal to the other 
boundaries (solid lines in Fig. ^,b). The boundaries 
of the model calculations were chosen to match those of 
electron micrographs of the actual SQUIDs measured. 
In choosing these boundary conditions we neglect the ef- 
fects of phase drops and supercurrent depression in the 
leads to and from the micro-SQUID. We argue that phase 
drops outside of the SQUID loop should have little effect 
on the critical current-flux characteristics which we are 
modelling, and that neglect of lead effects is therefore a 
good starting approximation. 



where i and j are indices labelling the 2-d matrices in the 
X and y directions respectively. Similarly, Eq. ^ becomes 

fi+l.J [fij + 2(/j+lj - fi,j)] + (y5j^lj/ij 
+ l [fiJ + 2(/i,j+l — fi,j)] + ipij-lfij 

+^,,, [-4/,,, - 2(/,+i,, - - 2(/,,,+i - = 0. (9) 

Using (in MKS units) | -0 m/2^oe^ X^f f where Ae// 
is the effective thin film penetration depth |l^ , the total 
current through the SQUID can be written as: 
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(10) 



c) Schemoi^c 

i 

s 

T 



d) Cross-sections 
1.0 



where Iq = dh/2^os^effj ^ thickness of the 

film from which the micro-SQUID is patterned. Super- 
current conservation requires that the sum over j in Eq. 

be independent of the value of i chosen. This was used 
as a self-consistency check of the solutions presented be- 
low. The difference equations Eq. are in the standard 
form 




FIG. 2. Calculations for the same geometry as Fig. Q but 
in this case for S/^=2.9, with a total phase drop ipa —4.4. In 
this case both the phase drop and superfluid density depres- 
sion extend throughout the SQUID structure. 
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Example solutions for and / using this model are 
shown in Fig. 's |l] and |[ Here the geometry is that of the 
Al micro-SQUID. The ipij and fij matrices both had 
102x110 elements, with 76 elements//im. The bridges 
were S'=289 nm long and 66 nm wide. The arms of the 
SQUID were 237 nm wide. Figure |^ shows the solution 
to Eq.'s H and ^ for S'/^=0.5, ips — tt. For a "short" 1-d 
micro-bridge {S/£^ << 1) with Lpg = tt, the solution for 
if has a step from to tt, and a depression in / to 0, 
with a characteristic width of ^, centered at the center 
of the bridge ||T^. The numerical solution for the full 
SQUID structure has a step in if and a depression in / 
near the center of the bridges, but there are significant 
changes well into the body of the SQUID, and / does not 
go completely to 0, meaning that Is does not go to zero, 
until is slightly above tt. These effects become more 
pronounced as 5/^ increases. Figure |2] shows the results 
for 5/^= 2.9, <y9s=4.4. This is the maximum value of 
for which a numerical solution to the Ginzburg-Landau 
equations could be found for S/£^ =2.9. In this case al- 
most half of the total phase drop occurs outside of the 
micro-bridge region. 
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FIG. 3. Calculated supercurrent-phase relationship for the 
Al SQUID, for values of S/^ appropriate for the data shown 
in Fig. I 



Supercurrent-phase relationships ls{}p) for this geom- 
etry are shown in Figure |[ This procedure is unable to 
find solutions numerically for the lower branch of the 
Is{^) characteristic when this characteristic is doubly 
valued [||,^. This will not affect our results, since we are 
comparing our modelling with the maximum supercur- 
rent at each value of applied field. This conclusion is sup- 
ported by our calculations of Ic{^) characteristics using 
the Likharev-Yakobson ^ Iciv') characteristics: only the 
upper branch is important for calculating the Ic{^) char- 
acteristics. The coherence lengths ^ in Fig. ^ were chosen 
to be appropriate for an Al micro-SQUID at the temper- 



atures for which the measurements of Fig. ^ were made: 
We use the dirty limit expression for the temperature de- 
pendent coherence length ^(T) = 0.855a/^oV(1 " T/T^) 
]l3[ , estimate the low temperature mean free path I in our 
films to be about lOnm from transport measurements, 
and take a value for of lOOnm and Tc—l.25K. 

A SQUID with non-sinusoidal current-phase {Ic{(p)) re- 
lationships will have critical current vs. flux interference 
patterns (/c($)) that are also non-standard. These can 
be modelled as follows. Assume that the SQUID has two 
arms labelled a and &, with total phase drops across the 
two arms ipa and cpb, inductances La and Lb, and micro- 
bridge supercurrents Isa = Lgi'Pa) and I^b = hdiVb)- 

The requirement of a single valued superconducting 
order parameter leads to 

27m ^ ip^- ipf,+ 2Tr4)e + PaQiVa) - f3bg{(pb), (12) 

where 0e is the externally applied flux <I>e divided by 
the superconducting flux quantum $o = h/2e, Pa ~ 
2TTLaIa/^o and f3b = 27rLf,/f,/$o- The total supercur- 
rent through the SQUID is 



(13) 



a) Experiment b) Short bridge c) Full model 
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FIG. 4. (a) Experimental measurements of the critical cur- 
rent of an Al micro-SQUID as a function of applied field, for 
selected temperatures, (b) Modelling in the short-bridge limit 
I=Iisiniy9, and (c) Modelling using the full Ginzburg-Landau 
calculations described in the text. 



The dependence of the critical current on applied field 
is determined most easily by assuming values for and 
one of the micro-bridge phase differences ifa or ipb, then 
varying the phase of the second micro-bridge until Eq. |l^ 
is satisfied. The values for ipa and ipb are then substituted 
into Eq. ^ to find the total current. The maximum value 
for / after repeating this procedure for all initial values 
of ifa and ipb is the critical current for that value of (/)e. 
The effect of the self-induced field |^] is expressed in Eq. 
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|T| via /3(r) = 27rL/c(T)/$o- /3(T) increases from in 
proportion to the critical current, since the inductance is 
purely geometric. In the case of the Al micro-SQUIDs, 
with a critical current of 170 /3(0.2 K) is 0.78; in 
the case of Nb, which have critical currents of 2000 /iA, 
(3(0.2 K) can be as high as 8.5. However, the differences 
in Ic are not enough to explain the differences in the 
Ic{^) characteristics between Al and Nb: they also arise 
from the shorter coherence length in Nb, which leads to 
more highly non-sinusoidal current-phase relationships in 
micro-SQUIDs made from Nb. 

Figure ^ shows critical current vs. applied magnetic 
field characteristics for our Al micro- SQUID. Figure ^ 
shows the prediction of Eq. |l^ and Eq. ^ assuming 
the standard Josephson current- phase relationship Ic=Ii 
sm{ip) for the micro-bridges. In this modelling the crit- 
ical currents of the bridges and the inductances of the 
arms of the SQUIDs were assumed to be symmetric and 
the critical currents were chosen to match the zero field 
critical current of the SQUID at each temperature. The 
inductances and Lf, were chosen to be 0.76pII each, 
half of the calculated total inductance L — 5/ioC/16, 
where C is the inner circumference of the SQUID |17||. 
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FIG. 5. Plot of the experimentally measured critical cur- 
rents at zero field (open circles) and at an applied flux of half 
the full modulation period (open diamonds), for an Al mi- 
cro-SQUID, as a function of temperature. The solid triangles 
are modelling as described in the text for the conventional 
sinusoidal current-phase relationship. The solid inverted tri- 
angles and squares are modelling using the Ginzburg-Landau 
calculations described in the text. 

Figure^ shows the predictions using the full GL calcu- 
lations for the current-phase relationship, with the same 
values for the inductances. For this modelling we take 
Icaiy^a) = Icai^b) = -^c(<y5)/2, whcrc Ic{(p)/h IS shown for 
the Al micro-SQUID in Fig. ^. In the determination of /q 
we assume \lff[T) = A2^^(0)/1 -t^,t = T/T^, d=38nm 
and use Ae//(()) as the sole fitting parameter. The best 
fit between experiment and modeUing for /c(T, $ = O) 



is for Ae//(0)=172 nm. This is to be compared with Al 
= 44 nm for bulk aluminum [T^ ]. The effective penetra- 
tion depth is often longer in thin films than in bulk, and 
can also be increased by impurity scattering ||l^,|l^ . Fur- 
ther, the effective thickness of the films will be reduced 
by oxidation. Once this scaling is done, the agreement 
between theory and experiment for /c(T, $ = 0) is good. 
The full model (Figure ^) fits the /c(T, $) experimental 
data well. This is made clear in the plot of /c(0), the 
critical current at zero applied field, and /c(<i>o/2), the 
critical current at the first minimum of the interference 
pattern, as a function of temperature, in Figure |5[ In this 
figure, the open circles are /c(0), and the open diamonds 
are /c($o/2). The solid triangles are the predictions for 
a sinusoidal current-phase relationship model as outlined 
above. The solid squares and inverted triangles are the 
predictions using the GL current-phase relationship de- 
scribed above. 
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FIG. 6. Calculated supercurrent-phase relationship for the 
Nb SQUID, for various values of S/^ appropriate for the data 
at selected temperatures shown in Fig. R 



Similar conclusions can be drawn from the data and 
modelling for a Nb micro- SQUID (Fig. I). In this case 
the /ij and (^i.j matrices were 100x100, with 67 ele- 
ments/micron. The micro-bridges were 184 nm long and 
100 nm wide, with the arms of the SQUID 285 nm wide. 
Iciv) characteristics for the values of S/^{T) appropri- 
ate for the data of Fig. |^ are shown in Fig. ^. In this 
case S/(,{T) is large, the variations of ip and / are spread 
throughout the SQUID structure, and /c(<y5) is nearly lin- 
ear, at all temperatures. For the calculations of the Ic{^) 
characteristics in Fig. the total inductance of the 
SQUID was taken to be 1.4 pH, the dirty limit expres- 
sion = 0.855(Co/)^/V(l - ty^^ was again used for 
the coherence length, with /=6.5nm ||l^, — SQnm [ p^ , 
Xlff{T) = A2^^(0)/l-^^ d=30nm, and = 8.23K. Fig- 
ure 0a shows experimental measurements of the depen- 
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dence of the SQUID critical current on applied magnetic 
field at selected temperatures. A multivalued critical cur- 
rent appears for currents higher than 0.4 and 0.5 mA for 
our given inductance. The sections of of the Ic{^) char- 
acteristics are nearly linear, with sharp discontinuities in 
the slopes at $ = (n + l/2)<i>o, n an integer. 

a) Experiment b) Short bridge c) Full model 

2.5 I • • • 1 I • • • 1 I • • • 1 
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FIG. 7. (a) Experimental measurements of the critical cur- 
rent of an Nb micro-SQUID as a function of applied field, for 
selected temperatures, (b) Modelling in the short-bridge limit 
I=Iisin<^, and (c) Modelling using the full Ginzburg-Landau 
calculations described in the text. 

Figure ^ shows modelling for an assumed sinusoidal 
current-phase relationship, using symmetric experimen- 
tally determined values for the micro-bridge critical cur- 
rents and symmetric inductances of 0.7pH for each arm of 
the SQUID. As for the Al micro-SQUID case, this mod- 
elling does poorly in describing the Ic{^) characteristics. 
Figure |^c shows modelling using the GL calculations ex- 
pressions for the micro-bridge current-phase relationship, 
with symmetric inductances of O.TpH in each arm, and a 
best fit value of Ae//(0)= 173 nm. This is to be compared 
with Ai,(0)= 44 nm for bulk Nb with a T^ of 9.26K 
As can be seen from Figures 0c and |[ the full modelling 
predicts the temperature dependence of the SQUID crit- 
ical current well, shows a triangular /c($) dependence, 
and does significantly better than that using a sinusoidal 
current-phase relationship for the modulation depth. 



DISCUSSION 

Our data shows that the modulation depth is reduced 
compared to the short bridge model as soon as the co- 
herence length becomes shorter than the bridge length. 
This effect is much more pronounced in the case of Nb as 
the intrinsic short coherence length of Nb leads rapidly 
to values of S/^ larger than 1. The model describes very 
well the overall lineshape in the case of both the Al micro- 



SQUIDs, as well as the Nb micro-SQUIDs, which show 
pronounced triangular /c($) characteristics. It is remark- 
able that the Nb micro-SQUIDs show quantum interfer- 
ence at all temperatures, even when the micro-bridges are 
much longer than the coherence length, and the phase 
gradients and supercurrent density depressions extend 
throughout the body of the SQUID. This implies that 
the heat dissipation that takes place when the SQUID 
enters the voltage state occurs well outside of the micro- 
bridge regions. 



Theory: 
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FIG. 8. Plot of the experimentally measured critical cur- 
rents at zero field (open circles) and at an appHed fiux of half 
the full modulation period (open diamonds), for a Nb mi- 
cro-SQUID, as a function of temperature. The solid triangles 
are modelling as described in the text for the conventional si- 
nusoidal current-phase relationship. The solid inverted trian- 
gles and squares are modelling using the full Ginzburg-Landau 
calculations for the current-phase relationship, as described in 
the text. 

How can the sensitivity of these micro-SQUIDs be im- 
proved? The Nb SQUIDs are limited in their sensitiv- 
ity by the short coherence length of Nb, as the modula- 
tion depth is maximal if the coherence length is equal to 
or longer than the bridge length. The coherence length 
could be increased by increasing the mean free path (epi- 
taxial films or thicker layers), but an increase in criti- 
cal current must be avoided. The performance of these 
micro-SQUIDs could be improved by significantly reduc- 
ing the length of the Dayem bridges, perhaps using novel 
scanning probe techniques . Decreasing thickness and 
width of the micro-bridge may lead to a diminution of 
the mean free path and thus to a diminution of the co- 
herence length. A further way to increase sensitivity of 
Al as well as Nb micro-SQUIDs may be the suppression 
of thermal hysteresis. A non-hysteretic behavior may be 
attained by fabricating the entire SQUID on a normal 
metal plane in order to remove the heat produced by the 
non local phase relaxation in the microSQUID. The non- 
hysteretic behavior will enhance the intrinsic bandwidth 
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and allow for a standard DC-SQUID detection scheme. 



CONCLUSION 

In summary, the temperature and field dependence of 
the critical currents in micro-SQUIDs can be understood 
by means of numerical calculations based on the phase 
dependence of the critical current predicted by a 2-d 
Ginzburg-Landau numerical calculation. It is important 
to do this calculation for the full micro-SQUID structure 
in many cases, because of the spreading of the variation 
in the superconducting phase and depression in the su- 
percurrent density beyond the micro-bridge region in this 
type of device. 

ACKNOWLEDGMENTS: 

J. K. would like to thank F. Tafuri for useful conversa- 
tions, and the Universite Joseph Fourier of Grenoble for 
their support during his visit in 1998, and A. Benoit is 
thanked for his encouragement of this work. 



[1] C. Chapelier, M. El Khatib, P. Perrier, A. Benoit 
and D. Mailly SQUID 91, Superconducting devices and 
Their Applications, edited by H. Koch and H. Liibbig 
(Springer- Verlag, Berlin, 1991) 286-291. 

[2] C. Chapelier, A. Benoit and D. Mailly, Phys. Rev. Lett., 
70 (1993) 2020. 

[3] W. Wernsdorfer, E. Bonet Orozco, K. Hasselbach, A. 
Benoit, B. Barbara, N. Demoncy, A. Loiseau, H. Pascard 
and D. Maillly, Phys. Rev. Lett., 70 (1997) 1791. 



[4; 

[5] 
i6l 

~] 

[8] 

[9: 

[10 

[11 

[12 
[13 
[14 

[15 

[16 
[17 
[18 
[19 



M. Jamet, W. Wernsdorfer, C. Thirion, D. Mailly, V. 
Dupuis, P. Melinon, and A. Perez, Phys. Rev. Lett. 86 
(2001) 4676. 

K. Hasselbach, C. Veauvy, D. Mailly, Physica C 332 
(2000) 140. 

C. Veauvy, K. Hasselbach, and D. Mailly, cond 



mat/01 10196 



T.A. Fulton, L. NN. Dunkleberger and R.C. Dynes,P/iys. 
Rev. B. 6, 855 (1972). 

K.K. Likharev and L.A. Yakobson, Sov. Phys. Tech. 
Phys., 20, 950 (1976) (Zh. Tekh. Fiz. 45, 1503-1509 
(1975). 

K.K. Likharev, Rev. Mod. Phys. 51, 101 (1979). 

E.P. Harris and R.B. Laibowitz, IEEE Trans. Magn. 

MAG-13, 724 (1977). 

W.-T. Tsang and T. Van Duzer, J. Appl. Phys. 47 (1976) 
2656. 

M. Faucher, T. Fournier, B. Pannetier, C. Thirion, W. 
Wernsdorfer, J.C. Villegier, and V. Bouchiat, Proceed- 
ings of SQUID 2001, to be published in Physica C. 
M. Tin- 
kham. Introduction to Superconductivity, (Mcgraw-Hill, 
New York, 1975). 

W.T. Vettering, S.A. Teukolsky, W.H. Press, and B.P. 
Flannery, Numerical Recipes Example Book [C], Second 
Edition, Cambridge, New York, 1988, p. 284. 
A. Baratoff, J. A. Blackburn, and B.B. Schwartz, Phys. 
Rev. Lett. 25 (1970) 1096. 

M.D. Maloney, F. de la Cruz, and M. Cardona, Phys. 
Rev. Lett. 5, 3558 (1972). 

Jacox and Ketchen, IEEE Trans. Magn. MAG-ll 
400. 

A. Barone and G. Paterno, Physics and Applications of 
the Josephson Effect, Wiley, New York (1982) p. 200. 
V. Bouchiat, M. Faucher, C. Thirion, W. Wernsdorfer, T. 
Fournier, and B. Pannetier, Appl. Phys. Lett. 79 (2001) 
123. 



7 



